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Abstract 

We calculate the Debye and Meissner masses of a gauge boson in a material consisting 
of two species of massless fermions that form a condensate of Cooper pairs. We perform 
the calculation as a function of temperature, for the cases of neutral Cooper pairs and 
charged Cooper pairs, and for a range of parameters including gapped quasiparticles, 
and ungapped quasiparticles with both quadratic and linear dispersion relations at low 
energy. 

Our results are relevant to the behavior of photons and gluons in the gapless color- 
flavor-locked phase of quark matter. We find that the photon's Meissner mass vanishes, 
and the Debye mass shows a non-monotonic temperature dependence, and at tempera- 
tures of order the pairing gap it drops to a minimum value of order ^Jly. times the quark 
chemical potential. We confirm previous claims that at zero temperature an imaginary 
Meissner mass can arise from a charged gapless condensate, and we find that at finite 
temperature this can also occur for a gapped condensate. 



1 Introduction 



In this paper we calculate the Debye and Meissner masses of a gauge boson propagating 
through a material consisting of two species of massless charged spin-^ fermions. We assume 
that, via an unspecified pointlike attractive interaction, these form an s-wave (rotationally 
invariant) condensate of Cooper pairs, which may be neutral or charged depending on the 
charges of the fermions. We allow the two species to have chemical potentials fl ± 6fi, and 
pairing gap parameter A, which is momentum-independent because the interaction is point- 
like. As one varies the chemical potential splitting 6fi of the two species, the spectrum of 
fermionic excitations (quasiparticles) changes dramatically: 

A> 6fi Gapped spectrum e.g.: Fig. [H dotted line 

A = 6fi Ungapped quadratic spectrum e.g.: Fig. [TJ solid line (1.1) 

A < 6fi Ungapped linear spectrum e.g.: Fig. [TJ dashed line 

We calculate the zero-momentum current- current correlation function (i.e. the Meissner and 
Debye masses of the corresponding gauge boson) for all these cases, for both a charged and a 
neutral condensate, as a function of temperature. 

To explain the motivation for this calculation, we give in section |21 a summary of the 
properties of ultra-dense quark matter, focussing in particular on the gapless color-fiavor- 
locked (gCFL) phase P, in which all these different cases occur^. Our calculation is described 
in sections El and |31 with technical details in an appendix. The results are presented and 
discussed in sections |S] and El with concluding remarks in section [7| 

2 Gapless color-flavor-locked quark matter 
2.1 Overview of the gCFL phase 

The behavior of matter at ultra-high density has been the subject of much theoretical work, 
and is gradually being constrained by experiments. It is generally agreed that at sufficiently 
high density, matter will be in a color-fiavor-locked (CFL) color-superconducting quark matter 
phase |2] (for reviews, see Ref. jS]). However, there are major questions about the next phase 
down in density. Recent work P suggests that when the density drops low enough so that the 
mass of the strange quark can no longer be neglected, there is a continuous phase transition 
from the CFL phase to a new gapless CFL (gCFL) phase. The two are very different: CFL 
quark matter is a transparent insulator, with no electrons jl]. Its only light mode is a neutral 
superfiuid mode associated with the breaking of baryon number. In contrast, gCFL quark 
matter is more like a metal: it has gapless quark modes, and electrons, as well as the superfiuid 
mode. 

There are many unanswered theoretical questions about gCFL quark matter. It may be 
modified by condensation of "kaons" j3] . Also, initial calculations indicate that some of the 
gluons have imaginary Meissner masses, indicating an instability towards the development of 
color currents [HI El El- Our results (section EJ confirm that imaginary Meissner masses are 
associated with gapless charged condensates. If the gCFL phase turns out to be stable after 

"'^The strong interaction that drives the quark pairing is not pointlikc. so the gCFL gap parameters are not 
momentum-independent. However, we do not expeet our results to be sensitive to this feature. 
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Figure 1: Dispersion relations for some of the quasiquarks in the gapless CFL phase. Note that 
there is a gapless mode with an approximately quadratic dispersion relation (solid line) as well as two 
gapless modes with more conventional linear dispersion relations (dashed line), and gapped modes. 
The labels give the color {r,g,b) and flavor {u,d,s) content of the quasiparticles. 

all, then an experimental question will arise: how could one detect gCFL quark matter in 
nature. The best opportunity for this is in the core of a neutron star, which achieves densities 
well above nuclear density, at low temperatures that permit color superconductivity. As we 
will discuss below, the most characteristic properties of gCFL quark matter are its transport 
properties. For example, the presence of a gCFL region will have a strong effect on the cooling 
of a neutron star ^j. Further progress will require the calculation of the interactions among 
the lightest excitations, and in this paper we will lay the groundwork for such investigations 
by performing a very general calculation of the behavior of zero-momentum gauge bosons 
coupled to both neutral and charged Cooper pair condensates. 

2.2 Quasiparticle dispersion relations in the gCFL phase 

In the gCFL phase, there is a condensate of Cooper pairs of quarks in the color-antisymmetric, 
flavor-antisymmetric, and Dirac-antisymmetric channel (there is also a insignificant color- 
symmetric flavor-symmetric component, which we neglect) 

(CCTsrf) ~ Aie-^ie,,i + A2e°^2^,fe2 + A3e°^3^,63 • (2.1) 

Here ip^ is a quark of color a = (r, g, b) and flavor a = {u, d, s). The gap parameters Ai, A2 
and A3 describe down-strange, up-strange and up-down Cooper pairs, respectively. 

The gCFL phase is named after its most striking characteristic: the presence of gapless 
modes in the spectrum of quark excitations above the color-superconducting ground state 
1)3.11) . These "gapless quasiquarks" were discussed in Ref. PP, and a sample dispersion relation 
plot is shown in Fig. ^ This was calculated using a Nambu-Jona-Lasinio model of the strong 
quark-quark interaction, with a quark chemical potential /2 = 500 MeV, strange quark mass 
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quasiparticles 


gapped 


gapped 


gapless 
(quadratic) 


gapless 
(linear) 



Table 1: The structure of the gCFL pairing pattern, which decomposes into one sector of 
three mutually paired species, and three sectors of two species that pair with each other. 
The Q charges of the quarks determine their electromagnetic interactions. Our two-species 
calculation can be applied to the gu-rd, rs-bu, and gs-bd sectors. 

rus = 200 MeV and pairing strength chosen such that at = the CFL gap would be 
Aq = 25 MeV pp. There are also gapped quasiquarks whose dispersion relations are not 
shown. 

The surprising feature of the gCFL phase is that there is a gapless mode with an approxi- 
mately quadratic dispersion relation E{p) oc {p — ppYi as "well as gapless modes with the more 
typical linear dispersion relation E{p) oc \p — Pf\^ and fully gapped modes. The quasiquark 
spectrum therefore includes all the cases listed in (jl.lj) . The presence of the quadratic gapless 
mode means that there is a lot of phase space at low energy, in fact the density of states 
for the quadratic mode diverges as E ^ We expect this to have dramatic effects on the 
transport properties, and it has already been pointed out that gCFL quark matter has an 
unusually large specific heat, which varies as cy oc \/T rather than the usual cy (x T that 
is associated with linear gapless dispersion relations. This will alter the late-time cooling of 
a neutron star in an observable way [Hj. To make further progress in studying the transport 
properties, it will be necessary to understand the dominant (electromagnetic) interactions of 
the lightest quasiparticles. For this we need the in-medium properties of the photon, which is 
why calculations of the Debye and Meissner mass of the photon are important. 

2.3 Relating our two-species calculation to the gCFL phase 

Our calculation explores the case of two species of fermions. This might seem to be a 
drastic simplification of the gCFL phase, in which there are nine species, but actually the 
gCFL pairing pattern decomposes into three sectors of two-species pairing and one of three- 
species pairing. Once pairing has occurred, the dominant interaction between the quasipar- 
ticles is mediated by the gauge boson of the remaining unbroken U{\)q gauge symmetry, 
whose gauge boson is mostly the original photon, with a small admixture of one of the glu- 
ons. The relevant quality of the quasiquark excitations is therefore their Q charge, given by 
Q = diag(|, — |, — |)flavor — diag(|, — |, — |)coior, not their electromagnetic charge. The pairing 
pattern and Q-charges of the quasiparticles are given in Table [U Each two-species sector has 
an average chemical potential /2 and a splitting Sfi that arises from the constraint of electrical 
neutrality and an approximate treatment of the strange quark mass, in which it is treated as 
a contribution — M^/(2/2) to the chemical potential for strangeness. This is known to be a 
good approximation PU} ITTj. 

From Table d we see that the transport properties will be dominated by the only elec- 
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tromagnetically interacting light degrees of freedom, namely the hu/rs quasiquarks (all the 
others are either gapped or Q-neutral) and the electrons. At any nonzero temperature there 
will be a nonzero density of these Q-charged particles which may lead to screening of the 
Q-electromagnetic fields. Because of their divergent density of states at low energy, we ex- 
pect the hu-rs quasiparticles to dominate this screening. The most direct application of the 
calculations in this paper is therefore to determine the in-medium properties of the photon 
in gCFL matter, taking into account the effect of the hu/rs quasiquarks. Since the hu and rs 
quarks have opposite charge, the relevant results are those for a neutral condensate, presented 
in sectional 

We should note, however, that our calculations of Debye and Meissner masses for charged 
condensates are also relevant: they can shed light on the behavior of the gauge bosons asso- 
ciated with the broken S'f/(3)coior generators, the gluons^. The diagonal color generators can 
be treated as U{1) gauge fields with their own charge assignments to the quarks, different 
from the Q charges. The fact that they are broken corresponds to the fact that some of the 
condensates have a net charge. Our two-species calculations show that, as one would expect, 
when the pairing species have charges that do not cancel each other there is a non-zero Meiss- 
ner mass for the gauge boson. When the pairing becomes gapless, we find that this Meissner 
mass becomes imaginary. This confirms existing zero-temperature calculations for the full 
two-flavor three-color (2SC/g2SC) and three-flavor three-color (CFL/gCFL) pairing patterns 
13 13 El- We find that even in the gapped case, turning on a temperature in the appropriate 
range can cause the Meissner mass to change from real to imaginary. 

3 Two-species pairing formalism 

In our analysis, we will treat two massless species of fermions (we will refer to them as 
"quarks") that pair to yield a condensate 

(^aC75^6) = Ae,b . (3.1) 

By comparison with Eq. (j2.H) and tabled we can see that the rd-gu, ru-bs, and gs-bd sectors 
of the CFL or gCFL quark condensate have this structure. 

Note that the Dirac charge conjugation matrix C does not connect left-handed with right- 
handed quarks: the pairing pattern is {'iPl^'l) and {ipB^n)- Similarly the gauge interactions 
preserve chiral symmetry, and do not couple left-handed quarks to right-handed quarks. A 
fermion mass term would couple left-handed to right-handed, but most existing treatments 
of the gCFL phase work to lowest order in the strange quark mass M^, including it via a 
chemical potential for strangeness = ~M^/P', which also does not couple left-handed 
quarks to right-handed quarks. To this level of approximation, then, we can treat the left- 
handed and right-handed quarks as completely decoupled from each other. We can therefore 
reduce the 4-dimensional Dirac space to a 2-dimensional Weyl space. To treat the quark- 
quark condensation, which violates fermion number and allows quarks to turn into antiquarks, 
we need to use Nambu-Gor'kov spinors, which incorporate particles and antiparticles into 
the same spinor, which doubles the size of our space. So we will work with 8-dimensional 

^The eighth gluon mixes with the photon, so the corresponding broken generator is also associated with a 
photon-ghion mixture. The other broken generators are entirely gluonic. 
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chiral spinors Xi which arise from a tensor product of the 2- dimensional Weyl space, the 
2-dimensional flavor space, and the Nambu-Gor'kov doubling. 
Ignoring electromagnetism, the Lagrangian for our quarks is 

^^ = \j xipY S-\p)x{p) dp , (3.2) 
where the inverse propagator is 

(p° + ;U) 1 ® 1 + p ■ (J ® 1 + 5/2 1 ® (Ts i A (72 ® CTi 



^ \^ -iA ® (Ti - /i) 1 ® 1 + p ■ ® 1 - 6^1 ® as 

(3.3) 

The Nambu-Gor'kov space is shown explicitly in the 2x2 structure of the matrix. In each 
entry, the first factor in the tensor product lives in the 2-dimensional Weyl (spin) space, and 
the second factor lives in the 2-dimensional flavor space. The only parameters are the average 
chemical potential /i, the chemical potential splitting 6n, and the pairing gap parameter A. 
The off-diagonal terms correspond to the quark pairing: they are proportional to A, with a 
Weyl factor of a2 because the Dirac matrix C = diag((T2, (T2) in the chiral basis, and a flavor 
factor of (Ti, indicating that the two flavors pair with each other, not with themselves. The 
on-diagonal terms are standard free fermion terms. 
The eight eigenvalues of are 

/ ± V(|^ ±^)2 + A2 ± 5fx. (3.4) 

The dispersion relations of the quasiquarks are given by the poles in the propagator, i.e. the 
zeros of det(S'~^). With the usual convention that negative energy states are filled, so that all 
excitations have positive energy, we find 



Eip)= V(|pl±/i)' + A2±5/i . (3.5) 

We see immediately how the relative sizes of A and determine the form of the quasiparticle 
spectrum, 

A > 6fi Gapped spectrum {rd/gu quarks in gCFL) 

A = 5ji Ungapped quadratic spectrum {hu/rs quarks in gCFL) (3.6) 

A < 5ji Ungapped linear spectrum [gs/hd quarks in gCFL) 

Actually, in gCFL the hu/rs is not precisely gapless, but it is very close: Sfi is slightly bigger 
than A, 

< A < A < /i , (3.7) 

and we will study how the Debye and Meissner masses vary as the temperature is scanned 
across the full range from T <^ 6fi — A to T ^ fi. Note, however, that because we do not 
include the physics that leads to the Cooper pairing, A is just a numerical parameter and 
it does not have the correct T-dependence that would send it to zero at ~ 0.57A. This 
means that our treatment of pairing is only valid for T <^ A (which is in any case the relevant 
regime for quark matter in neutron stars) so insofar as our results depend on A, they are only 
valid for T < A. 



4 Gauge boson self-energy 



We now discuss the in-medium properties of a gauge boson. These depend on the charges of 
the two quark species. In general their charges could be {qi,q2) = Q ^ ^1, and it turns out 
that these contributions decouple from one another, so that for both the Debye and Meissner 
masses we find 

= g'M,L,g,d + (5g)'Mncutral (4-1) 

Where Mcharged is the result for a charged condensate in which both flavors have the same 
charge, and Mneutrai is the result for a neutral condensate in which the two flavors have 
opposite charges. Because of this decoupling, we will consider separately the case of a neutral 
condensate with (gi, ^2) = (1, —1) and a charged condensate with (gi, 52) = (1, !)• 

In the Nambu-Gor'kov formalism, the covariant coupling of the fermions to the gauge 
boson takes the form 

^ex^r'^X (4.2) 
where the gauge coupling is e and depends on the charges of the fermions: 

pM ^ ( crf^^a^ 

neutral Q -a'''^ ^ (T3 

PM ^ [ a^'®! \ (4.3) 

charged g -(T^^ <S> 1 J ' 

= (l,a). 

To lowest order in the gauge coupling, the gauge boson self-energy at external momentum q 
is 

KIM) = e'J -0^jT[T^S{k)T'^S{k-q)]. (4.4) 

This corresponds to the Feynman diagram shown in Fig. |2l It is to be evaluated using the 
fermion propagators S obtained by inverting (j3.3j) . The technical details are given in the 
appendix. Note that our expression ()4.4|1 lacks a factor of | when compared with Eq. (20) of 
Ref. jT2j or Eq. (35) of Ref. [B]. This is because we only include one chirality of fermion in 
our formalism, so we multiply our result by 2 to obtain the value of Ilbare for a Dirac fermion 
with both chiralities. 

For the transport properties of gCFL quark matter, we are interested in low momentum 
properties, corresponding to the limit g — 0. For applications to neutron stars we are inter- 
ested in low temperatures T ^ 0. We also expect singular behavior as the quark spectrum 
goes quadratically gapless {6fi A) as seen in the bu/rs modes in gCFL. We have to be 
careful about the ordering of all these limits. We will take g ^ first, and then study a range 
of T, both above and below \6fi — A\. Our calculations of Ilbare will yield divergent integrals 
such as (IA.17j) . To obtain physical results we require a regularization and renormalization 
prescription. We use a momentum cutoff, and subtract off the vacuum contribution, so that 
our renormalized results contain only the parts that depend on /2, 5/1, and T, 

n(T, /i, A) = nbare(T, /i, A) - nbarc(0, 0, 0, A) . (4.5) 

The renormalization subtraction term uses the same gap parameter A as the bare contribution. 
Ultimately, the justification for this is that we then get the correct value for the Meissner mass 
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Figure 2: The Feynman diagram for the lowest-order contribution to the photon self-energy. The 
internal lines are fermion propagators, obtained by inverting l\3.'d\i . 

in the presence of a neutral condensate, namely zero. Other calculations in the literature 
O El use a different prescription, n(T, /i, 5/1, A) = nbarclT", /i, 5/1, A) — nbarc(0, 0, 0, 0), which 
in our calculation would yield a spurious contribution of order e^A^ to the Meissner mass. 
Such a contribution can be seen in Rischke's result for the A3 gluon mass in two-flavor quark 
matter (Eq. (115) of Ref. [12] )• Rischke guessed that this contribution would be cancelled 
by other contributions that had been neglected in his calculation. However we do not make 
any approximations (unlike Rischke we do not try to combine a realistic pairing mechanism 
with our gauge boson mass calculation, so our computations are simpler than his) and we can 
see that this term is not cancelled. We think that our prescription is the correct one for the 
situation that we study, but we do not claim to have provided an a jjnon justification for it. 

5 Debye and Meissner mass for a neutral condensate 
5.1 Results 

The Debye and Meissner masses are defined by the static long-distance limit of the self energy 
n'^'^(g). (Analysis of the self energy for non-zero q yields interesting information about how 
the photons resolve the Cooper pairs ^|, but we do not attempt such an analysis here.) 
Taking g ^ in equations ()A.15|) and ()A.17|) and using (j4.5|) . we obtain 

limn°°(g) = ^2fdkei^^[n^{E^)+n.{E^)] + ^ [n^{E_) + n^E^)]^ , (5.1) 
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and 
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where 



B±(*-) = ^{k ± lip + A2 , 
".W - (e.p(^) + l) 



-1 



(5.2) 

(5.3) 
(5.4) 



The integral ()5.2p evaluates to zero as one would expect: a neutral condensate does not 
break gauge symmetries, so the Meissner mass is zero, 



M^ = ^lim(5,,-g,g,)ff^(g) = 0. 



(5.5) 



Note that in obtaining this result it was crucial that we used the correct renormalization 
subtraction (j4.5|) . 

The integral in ()5.ip can be evaluated numerically, and also analytically in certain limits. 
For the Debye mass Md, defined by M|, = — limg^o n°°(g), we find 



2 1 



2 1 



6fi 



271 A A-^M 

^ e I' . 
T 



where the numerical constant C is 

c = 



A < T, /i (a) 

A| <T< A,5^<^ (b) 

, T < A|, A < 5/i < /i (c) 

T < A|, < A < ^ (d) 



(5.6) 



dx- 



oo (e-' + 1)' 



0.673718 . 



(5.7) 



In the zero temperature limit this becomes 

2e^ 



MjjiT = 0) 



(/x^ + 5/i^-A^ 



5ii 



^((5/i - A) 



(5.8) 



Q 




Figure 3: How the dispersion relations of the quasiparticles look at energies of the same order as the 
temperature, for temperatures coresponding to the cases distinguished in Eq. (|5.6() . The temperature 
decreases from (a) to (b) to (c)/(d), and the energy scale stretches accordingly. In each case the 
shaded region indicates energies E < T. 



5.2 Discussion 



Our result for the Debye mass of a photon passing through a neutral condensate of charged 
quarks shows a subtle interplay of the limits of small T and small \6fi — A\. To understand 
why the various limits behave so differently, it is useful to recall how the dispersion relations 
look at energy of order T for the different ranges of temperature distinguished in Eq. fl5.6|) . 
These are shown in Fig. El We will discuss each of the four regimes, bearing in mind that the 
physical application of our result is to a photon in gCFL quark matter, in which the photon's 
in-medium properties are dominated by the bu/rs quarks with their near-quadratic gapless 
dispersion relation. 

Case (a) corresponds to high temperatures, where the structure in the dispersion relations 
at scales of order A, 5/i is invisible, so the system behaves like free particles with chemical 
potential p.. As noted in section HI our treatment of the pairing is not valid in this range, but 
the result turns out to be independent of pairing (i.e. of A). Because we have two species 
we get double the standard one-species result, which is M|, = e^(/i^/7r^ + (Ref. [13] 

Eq. (6.103), in which = 2m'^, see Eq. (7.125)). 
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Figure 4: The Debye mass as a function of temperature, for p, = 500 MeV, 5fi ^ A = 25 MeV, 
(5/U — A = 10"^ MeV. Tlie lieavy dashed line is the numerical result (|5.1j) . The solid lines show the 
analytic approximations in various regimes of Eq. 1)5. 6|) . (a), (b), (c). The temperature of a neutron 
star varies from ~ 10 MeV in the first seconds after the supernova to ~ eV after hundreds of millions 
of years. 

Case (b) corresponds to temperatures much less than 6fi and A, but much greater than their 
difference. This is the typical situation for the bu/rs quasiparticles of gCFL matter in a 
neutron star that is not very old. At energies of the order of T the dispersion relation looks 
as if it barely touches the momentum axis at a double zero (compare the bu/rs line in Fig. P). 
This gives the unusual behavior M|, oc I/VT: the Debye mass increases as T decreases. 

Case (c) is for temperatures far below the splitting, in the case where 6fi > A. This corre- 
sponds to the bu/rs quasiparticles of gCFL quark matter in a very old, cold neutron star. Now 
a typical thermal fluctuation can resolve the apparent double zero into two separate zeroes. 
The Debye mass levels out at a constant value as it would for free particles. If A ^ ^ /i 
then this constant value is the same as in case (a). 

Case (d) is at a similar temperature to case (c), but for a system where 5^ < A, so the 
dispersion relation never actually drops to zero. Relative to a typical thermal fluctuation the 
quasiparticle gap is a large energy barrier, and Mo drops to zero very rapidly with decreasing 
T. 

This gives us the temperature dependence for M^) shown in Fig. |3 We have chosen values 
for the parameters that are appropriate for gCFL matter in a neutron star: e^/(47r) = 1/137, 
/i = 500 MeV, ^ A = 25 MeV, 5/x - A = 10"^ MeV [I]. We evaluate the integral 
numerically, and compare with the analytic approximations of ()5.()|) . Note that Mo reaches 
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Figure 5: The square of the Meissner mass mass as a function of the gap parameter A at various 
temperatures, for p, = 5, (5/i = 1 in arbitrary units. Note how the zero-temperature discontinuity is 
smoothed out at T > 0. This means that in the gapped case (A > turning on a temperature in 
the appropriate range can cause the the Meissner mass to become imaginary. 



a very large value at T = 0. This is because at T = 0, Md diverges as 6j2/ \/ — when 



6 Debye and Meissner mass for a charged condensate 
6.1 Results 

From ()A.21j) and (IA.22jl . using (jUSJ, we get 
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Figure 6: The square of the Meissner mass mass as a function of temperature for various values of 
the gap parameter, with p, = 5, Sfi = 1 m arbitrary units. As in Fig. I^J we find that for the gapped 
case (A > 5//), M^^ can go from a positive value at T = to a negative value when T > A — 6fi. The 
curves stop when T reaches 0.5A, at which point the temperature-dependence of A can no longer 
be neglected. 



At arbitrary temperature these integrals can be evaluated numerically. In the zero temperature 
limit they reduce to 

MUT^O) = g^jf.^ + 0(6^ - A) + A' log C'" ' 'Z^"' ' ^ )) } 

(6.3) 

When 5/i < A we recover the well-known result for gluons in the CFL phase, MIj/M"^ = 1/3 

unmni. 

6.2 Discussion 

For a given we can calculate the Debye and Meissner masses as a function of temperature 
T and the gap parameter A. The Debye mass is large, Mr, > /i ^ A, so we do not discuss 
it in detail. In the case of the Meissner mass, we can see that at T = we get an imaginary 
value when > A, i.e. when the system is gapless. This is reminiscent of results found by 
Shovkovy and Huang and others for gluons in the more complicated cases of g2SC and gCFL 
quark matter jHl 13 E] • The square of the Meissner mass as a function of temperature is shown 
in Figs. El and ini To make these plots we set fl = 5, 6fi = 1, in arbitrary energy units, and 
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used e^/(47r) = 1/137 as in section |3| The solid "T = 0" curve in Fig. El shows the result of 
Eq. (j6.3|) : the Meissner mass Mm is zero when there is no pairing (A = 0) ^ . When A < 
MIj is negative, and diverges to —oo as A ^ At A = 5/1, jumps discontinuously to 
the positive value le^Jj^j (Svr^) that is characteristic of the spontaneous breaking of a gauged 
symmetry by a charged condensate. It is clear that the limit A — 5/1 is singular, with very 
different behavior according to whether it is taken from above or below. This is quite natural, 
since for A > 6fi the system is always gapped, whereas for A ^ 6fi it is gapless. We expect 
that this singular behavior will be smoothed out at any finite nonzero temperature, and this 
is indeed the case. The curves for T = 0.1 and T = 0.3 in Fig. El show a smooth transition, 
over a range 6 A ~ T, from the gapless behavior to the gapped behavior. 

Essentially the same information is presented in a different way in Fig. IHl where we show 
the dependence on T for a range of different values of A. Again, we fix /i = 5, Sfi = 1. 
For A > Sfi (gapped system) we see Mfj tends to a positive constant as T — > 0, whereas 
for A < 6n we find that tends to a negative value that is large for A just below 5fi, 
but tends to zero as A — > 0. It is also apparent that the smoothing effect of the non-zero 
temperature can change a positive into a negative one. This effect is also visible in Fig.El 
and is quite reasonable, given the picture outlined in Fig. |21 If A > 6fi then the dispersion 
relations are gapped, with positive at T = 0; but if T > \6fi — A\ then to within 
the natural energy resolution (T, approximately) the dispersion relation looks quadratically 
gapless and the Meissner mass will become imaginary. It should be remembered that we 
have not included any of the pairing dynamics in this calculation, so our gap parameters are 
independent of temperature. The curves in Figs. Eland El stop when T reaches 0.5 A, at which 
point the temperature-dependence of A can no longer be neglected. 

7 Conclusions 

We have calculated the Debye and Meissner masses of a gauge boson in the presence of a 
condensate of Cooper pairs that involve two species of massless charged spin-| fermions. We 
did not specify any particular pairing mechanism, but simply parameterized the dispersion 
relations of the quasiparticles using a momentum-independent gap parameter A, and individ- 
ual chemical potentials /i ± 5/i for the two species. We allowed the fermions to have arbitrary 
charge in their coupling to the gauge boson, but found that this reduced to two elementary 
neutral condensate (charges (-1-1, —1)), and a charged condensate (charges (-1-1, +1)). 
Our results for the neutral condensate are presented in sectionEland our results for the charged 
condensate in section El 

For the neutral condensate, our results give the in-medium behavior of a very low energy 
photon (technically, it is the massless Q gauge boson that is predominantly the photon with 
a small admixture of gluon) in the gCFL phase of quark matter. This is dominated by the 
gapless charged excitations, the hu/rs quasiquarks, which form a two-species system of the 
type that we studied, where the strange quark mass has been treated in lowest order as a 
contribution to the chemical potential for strangeness. Because the photon coupling is weak, 
our calculation, which includes only the leading order diagram (Fig. |2I), gives the dominant 

^Technically our renormalization condition (|4.5(l corresponds to throwing away the part of the A- 
dependence that is independent of /i, T, Sji. But this is always zero because free fermions at T = /i = 
have no Meissner mass. So it is legitimate to calculate the A-dependence in our framework. 
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contribution. We find that the Debye mass shows a non-monotonic temperature dependence, 
dropping to a minimum value of order -y/a/i at temperatures somewhat below the pairing 
gap. We find that the Meissner mass is zero, as expected for a condensate that does not 
break the gauge symmetry, and we note that to obtain this result it was necessary to use the 
renormalization subtraction ()4.5|) . which differs from the one used in the existing literature 
on quark matter. 

For the charged condensate, our results give some insight into the in-medium behavior of 
gluons in color-superconducting phases. The color- diagonal gluons can be treated as photons 
with appropriate charge assignments to the various quarks so that the condensates in all three 
2x2 sectors of the color- flavor- locked pairing pattern (Table ^) carry net color charge, except 
that the rd-gu condensate is neutral to the A3 gluon. We find that there is always a large 
Debye mass, of order /i^ or greater. At zero temperature, the square of the Meissner mass is 
negative for < A ^ Sfi, diverging to —00 as A — > Sfi~ , but then jumps discontinuously to a 
positive value for A > 6fi. The Meissner mass is therefore imaginary whenever the quasiquark 
spectrum is gapless, and positive when the spectrum is gapped. At T > this discontinuity 
is smoothed out, with interesting consequences. For gapped systems (A > 6fj,) the Meissner 
mass is real at T = 0, but it becomes imaginary when T > A — as the temperature-smeared 
dispersion relation is then indistinguishable from a gapless one. 

Our charged-condensate results confirm the essential conclusion of Refs. [B] and [7j, that 
charged condensates with gapless excitations are associated with imaginary Meissner mass. In 
particular, our results for the Meissner mass agree well with those obtained for the diagonal 
"8 gluon" (the combination of a gluon and a photon that is orthogonal to Q) in Ref. [0]. It 
is interesting to note that Ref. [Oj also finds imaginary Meissner masses for the off-diagonal 
gluons in the g2SC phase, even in situations where no quasiquark modes are gapless. We 
cannot offer any insight into this because we only treat diagonal gauge bosons^. Also, it is 
curious that Ref. [7], which treats the gCFL phase, only finds an imaginary Meissner mass 
for the off-diagonal (Ai and A2) gluons. Although our 2-species calculation is not directly 
applicable to the gCFL case, which has multiple pairing sectors (table [T]) coupled together 
by neutrality constraints, our results would lead us to expect that they should have found 
imaginary Meissner masses for the diagonal A3 and 8 gluons, which both couple to the gapless 
bu-rs and gs-bd condensates. 

Although we have couched our discussion in terms of gauge boson masses, the quantity 
that we calculated, the low-momentum current- current two-point function, also has physical 
meaning if the currents in question are not coupled to gauge fields. In this case M|^/(e^A^) 
is the coefficient of the gradient term in the effective theory of small fluctuations around the 
ground-state condensate. The fact that we find a negative value when the quasiparticles are 
gapless indicates an instability towards spontaneous breaking of translational invariance. The 
nature of the true ground state in gapless color superconductors remains unknown: it could 
be a mixed phase ^Hl or a crystalline (LOFF) phase jH]. Since the two-species system shows 
this instability, it may provide a convenient toy model for investigating the nature of the true 
ground state. 

^Our formalism can easily accomodate an off-diagonal coupling. We have performed preliminary calcula- 
tions for this case, but find an unphysical 5/x-dependcnt logarithmic divergence. This was present but discarded 
in the recent 2-flavor (g2SC) calculation and is also thought to occur in the analogous gCFL calculation 
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Appendix A Calculational details 
A.l Fermion propagator 

We calculate the fermion propagator by inverting in ()3.3|) . To doing so, let us write S as 
a 2 X 2 matrix 

Since S is Hermitian, we have 

^12 = Si. (A.2) 
We may expand each matrix element in 5* as 

Su = -iA[a{p)li^l + b{p)lCSa3 + c{p)p-a®l + d{p)p-a^as]{cT2^ai); 

S22 = a{p)l ^ 1 + (3{p)l ^ as + -f{p)p ■ a'^ ^ 1 + Vip)p ■ ® (^3- (A.3) 

After a lengthy calculation, we have 

1 1 1 1 
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(A.6) 



where 

Di{p) = {p' + 5fiY-{\p\+fif-A' 

D,{p) = {p^ + 6f,Y-{\p\-fir-A' 

Dsip) = (/-5/i)^-(b1+/i)^-A^ 

D,{p) = (p^-Sf^f-{\p\-fif-A\ (A.7) 



A. 2 Matsubara frequency sums 



To compute n(^^j.^,, we need do evaluate the 4-dimensional momentum-space integral At finite 
temperature, the ko integral becomes a Matsubara frequency sum, in which we replace /c° by 
iun, and J dk^ by 27rT^^ . Because quarks are fermions, they obey antiperiodic temporal 
boundary conditions, and the Matsubara frequencies are 



= — (2n + l)7r, neZ. 

To evaluate the frequency sums, we use the following result (Eq. (5.77) in Ref. 

1 1 



1 



AE1E2 
+ 



f_{E,)-f4E2) U{E,)-U{E,) 



iuj + El — E2 

i-/4^i)-/+(^2) 



iu — El + E2 

I- U{Ei)- f^{E2) 



iu + Ei + E2 iuJ — El 

where with positive E, the functions f± are defined as 
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From Eq. (5.76) in Ref. 

1 



we obtain another useful formula 

1 f^SiEi) - f4s2E2) 



(A.8) 



iuJn + + SiEi i{uJn — Cij) + 5/i + 52-^2 + SiEi — S2E2 



(A.9) 
(A.IO) 

Si,S2 = ±1. 

(A.ll) 
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To apply these formulas directly, we need to use the method of partial fractions to convert 
expressions with fc^ in the numerator to expressions with fc" appearing only in the denominator. 
A typical example is 
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(A.12) 



With the integrand in this form, the frequency sum can be evaluated. 
A. 3 Neutral condensate 

We substitute the fermion propagator into the expression for the photon self-energy ()4.4|1 with 
^neutral (j4.H|l . The zcro-zcro component is 
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where — >■ ±Sn){fi — ±/i) means terms as the first two terms with 6fi replaced by ±5/i 
and fl replaced by ±/i. 

The spatial components of the photon self-energy are 
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where (/i ^ — /i) and — * — means a duplication of all previous terms with p, replaced 
by —fi or 6fi replaced by — 

Performing the frequency sums using ()A.9|) and ()A.11|) . we obtain the zero- zero component 
of the photon self-energy: 
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For the spatial components, we obtain 
^2 r 



e 
~4 



1 5*^' [l - L (k-q) + k\k^qy + {k^qyy 



ijl 



k 



ik^q)i] ] 



X 
X 
X 
X 



l + f^S{k,q) u^Ak,q)+ l-r:!+{k,q) v+,+ {k,q) 



d2) 



l-k-{k-q) + k\k-qy + {k-qyy 



+ ie 



ijl 



h - {k-q)i 



l + k-{k-q) - k\k-qy + {k-qyy 



ijl 



h + {k-q)i 



l + ftl{k,q) u+Ak.q)+ l-ftUk,q) v+^k^q) 



c{2) 



l + k-{k-q) - k\k-qy + {k-qyy 



ijl 



ki + {k-q)i 



f^^likq) v.Ak,q) 



where 



fSik,q) 



r\k\ + fi) { s\k — q\ + fi) — 



(A.16) 



(A. 17) 

(A.18) 
for the 



Erik)E,ik-q) 

Now we can send the external momentum g ^ 0, to obtain integral expressions 
Debye and Meissner masses (see section E}. 

A. 4 Charged condensate 

In the formalism we are using here, the only difference between the neutral case and charged 
case is the sign in front of A^ in the numerators. If we change all the signs of A^ in the 
numerators in the previous subsection and leave the denominators unchanged, we will obtain 
the correct formulas for the photon self energy in the presence of a charged condensate. 
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The zero- zero component of the photon self-energy is 
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The spatial components of the photon self-energy are 
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To perform the frequency sum we can take the neutral-condensate results and interchange 
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Now we can send the external momentum g ^ 0, to obtain integral expressions for the Debye 
and Meissner masses (see section IH)). 
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